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We present measurements of lepton {£) angular distributions in tt — > W + bW~b — > V~vbtTvb de- 
cays produced in pp collisions at a center-of-mass energy of yfs = 1.96 TeV, where I is an elec- 
tron or muon. Using data corresponding to an integrated luminosity of 5.4 fb _1 , collected with 
the DO detector at the Fermilab Tevatron Collider, we find that the angular distributions of t~ 
relative to anti-protons and l + relative to protons are in agreement with each other. Combin- 
ing the two distributions and correcting for detector acceptance we obtain the forward-backward 
asymmetry A FB = (5.8 ± 5.l(stat) ± 1.3(syst))%, compared to the standard model prediction of 
A FB (predicted) = (4.7±0.1)%. This result is further combined with the measurement based on the 
analysis of the £+jets final state to obtain Ap B = (11.8 ± 3.2)%. Furthermore, we present a first 
study of the top-quark polarization. 

PACS numbers: 14.65.Ha, 12.38.Qk, 13.85.Qk, 11.30.Er 



To check the validity of the standard model (SM) of 
elementary particle physics and to search for possible ex- 
tensions, we measure the properties of the top (t) quark. 
At leading order (LO) in perturbative quantum chromo- 
dynamics (QCD), production of pairs through quark- 
antiquark (qq) annihilation is expected to be forward- 
backward (FB) symmetric in the center-of-mass frame. 
At next-to-leading order (NLO) QCD, interference leads 
to a positive FB asymmetry, which implies that the top 
(antitop) quark is emitted with higher probability in the 
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direction of the incoming quark (antiquark). Top pair 
production through gluon-gluon fusion does not lead to 
such asymmetry 

SM predictions for the FB asymmetry can be modified 
by processes beyond the SM [1 J2| , such as contributions 
from hypothesized axigluons 3|], Z' or W bosons [H, 
and new scalars These sources of physics beyond 
the SM modify observables sensitive to the top quark 
polarization [6|. 

At the Tevatron pp collider, with y/s = 1.96 TeV, tt 
production is dominated by qq annihilation. The direc- 
tion of the incoming quark (antiquark) therefore coin- 
cides with the direction of the proton (antiproton). The 
t quark (antiquark) is more likely to be emitted in the 
direction of the incoming quark (antiquark) than in the 
opposite direction. This FB asymmetry in tt production 
can also be observed through the t and i decay prod- 
ucts, for example, in the distributions of charged leptons 
(I = e, n) from t -> W+b -> i+vb and t -t W-~b -> i~vb 
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decays. 

The CDF and DO Collaborations have previously per- 
formed measurements of the FB asymmetry in ti decays 
to £ + jcts final states containing exactly one leptomiets 
and an imbalance in transverse momentum (Ipr) |2ri2|- 
The asymmetries reported by both Tevatron collabora- 
tions are larger than predicted in NLO QCD. The asym- 
metry in CDF data at large values of ti invariant mass 
(m t > 450 GcV) differs by more than three standard de- 
viations (SD) from the NLO prediction Q. The DO data 
show no significant excess in this mass range. Defining 
a ti asymmetry based on the pseudorapidity, r\ (loj . of 
the charged lepton, DO finds a significant deviation from 
NLO QCD predictions of the order of three SD [§]. The 
ATLAS and CMS collaborations have performed mea- 
surements of the difference in angular distributions be- 
tween top quarks and antiquarks in the €+jets final state 
using asymmetries based on the top quark and antiquark 
rapidities l_l|, LL2J and pseudorapidities [12J . The results 
are consistent with the SM expectations. 

In this Letter, we present six measurements of lcptonic 
FB asymmetries in pp collisions at y/s = 1.96 TcV, us- 
ing data corresponding to an integrated luminosity of 
5.4 fb _1 , collected with the DO detector in Run II of 
the Fermilab Tevatron Collider. We use ti candidates 
in dilepton final states, where the W bosons from t and i 
decays both decay into ev e , fiv^, or tv t , and the r lepton 
decays leptonically (r — > lve.v T ). We calculate asymme- 
tries based on the pseudorapidity and charge of the elec- 
trons or muons. These asymmetries are determined from 
the angles of the charged leptons, which are measured 
with high resolution. These measurements have the ad- 
vantage that a full reconstruction of the ti event is not 
required. In addition, we combine this measurement of 
the FB asymmetry with the DO measurement performed 
in £ + jets final states @. Furthermore, we present a first 
study of the longitudinal polarization of the top quark. 

A description of the DO detector can be found in fl3j |. 
The selection criteria and object identification of the 
dilepton (ee, en, nu) decay channels follow those de- 
scribed in Ref. )14| . To enrich the sample in ti events, 
we require two isolated, oppositely charged leptons with 
transverse momentum pt > 15 GeV and at least two 
jets with pt > 20 GeV and detector pseudorapidity 
|?7dot| < 2.5 fioj . For the ep, channel we require that 
Ht (defined as the scalar sum of the larger of the two 
lcpton-pT values and the scalar pt of each of the two 
most energetic jets) be greater than 110 GeV. For ee and 
fifi events we compute a likelihood for the significance of 
JpT [HI ) based on the probability distribution calculated 
from the value of IpT and the lepton and jet energy res- 
olutions. We require this likelihood to exceed the value 
typical for background events. We find that only the 
fi/j, channel benefits from an additional restriction on I^t 
and, to increase signal purity, we therefore require Ifir 
> 40 GeV for the p,p, final state. We select a ti sample 



with a signal to background ratio of 3.2, 3.7 and 0.9 in 
the ee, e/x and fifj, final states, respectively. 

To simulate ti production, the MC@NLO [l6| generator 
is used assuming m t = 172.5 GcV. The production of top 
quarks is simulated at NLO, while the decay is simulated 
only at LO. To include full NLO QCD corrections to both 
production and decay as well as mixed QCD and quan- 
tum electrodynamic corrections and mixed QCD and 
weak corrections to the production amplitudes (denoted 
by "QCD+EW"), we simultaneously correct the nor- 
malized lepton and antilepton rapidity distributions in 
MC@NLO using the predictions of Ref. [I3| • HERWIG [I3| 
is used to simulate fragmentation, hadronization and de- 
cays of short-lived particles, and the generated events 
are processed through a full detector simulation using 
GEANT [l]|. The Monte Carlo (MC) events arc overlaid 
with data from random bunch crossings to model the ef- 
fect of detector noise and additional pp interactions. The 
same reconstruction programs are then applied to data 
and MC events. The background in the dilepton channel 
arises from Z/j* — > l + (r and diboson events (WW, WZ 
and ZZ) with associated jets, from instrumental back- 
ground where a jet is misidentified as a lepton, and from 
heavy quarks that decay into leptons that pass isolation 
requirements. A detailed description of these processes 
and their generation can be found in Ref. (20| . 

Leptons are reconstructed with excellent resolution on 
the measurements of their angles and electric charge. 
In contrast, it is challenging to reconstruct the four- 
momenta of the t and i quarks, since the kinematics is 
undcrconstrained because of the two neutrinos in the fi- 
nal state. Rather than reconstructing the t and i four- 
momenta, as in Refs. we measure observables cor- 
related to the FB asymmetry, which depend solely on 
the r\ and electric charge of the lepton £, as proposed in 
Ref. Q. The asymmetry for leptons is defined as: 
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N t+ (r)>0)-N t -(ri>0) 
N i+ (rt>0)+N t -(v>oy 



(1) 



where N#- (rj) and Ng+ (rj) correspond to the number of 
leptons and antilcptons as a function of r/, respectively. 
If CP invariance holds in ti production and decay, then 
Ne+(r]) = Nf-(—rf), and Ap B defines the FB asymmetry 
for both leptons: 



1 FB 



N e ±{7 1 >0)-N e ±(r ] <0) 
N i ±{7 1 >0) + N e ±(r,<0)' 



(2) 



The asymmetries Ap B and Ap B are statistically indepen- 
dent and opposite. We can therefore combine the asym- 
metries for £ + and l~ by multiplying r\ with the charge 
Q of each lepton: 



A 



FB 



N t {Q ■ r) > 0) - NjjQ ■ y < 0) 
N e {Q ■ rj > 0) + N e (Q ■ V < oy 



(3) 
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FIG. 1: Pseudorapidity distributions of the charged leptons 
for the combination of the ee, e/j and /j/j final states after 
the selection criteria have been applied. The r\ distribution 
of positively (a) and negatively (b) charged leptons, the dis- 
tribution of Q ■ 77 (c) and the distribution of A?; = ry f+ — r\ t - 
(d) are shown. The vertical error bars indicate the statistical 
uncertainty. The ti contribution is normalized to the data 
after background subtraction. 



In analogy to the FB asymmetry for t and t quarks, 
we define an angular asymmetry for leptons: 
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N(Ar) > 0) - N{Arj < 0) 
N(Ar] > 0)+N{Ar] < 0) ' 



(4) 



where A77 = rji+ — . The asymmetry corresponds 
to a longitudinal asymmetry in spin orientation relative 
to the proton beam direction. It is defined as: 
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CP 



N i+ {r 1 >0)-N i -(r ] <0) 
N t+ (ri>0) + Nt-{n<0)' 



(5) 



This asymmetry is sensitive to s-channel exchanges of 
heavy non-scalar resonances with CP- violating couplings 
to quarks, but not to possible P and CP-violating effects 
from an s-channel exchange of Higgs bosons . 

The asymmetries are measured in four ways using 77 
and Q of the leptons: separate 77 distributions for (i) £ + 
and (ii) £~ , (hi) the charge-signed pseudorapidity, Q ■ 77, 
and (iv) A77. They are presented in Fig. [I] To extract the 
asymmetries for ti events from the distributions shown in 
Fig. [U we subtract the background and then correct for 
effects from event reconstruction and acceptance. The 
correction for detector acceptance is performed by multi- 
plying the background-subtracted number of events with 
the inverse of the selection efficiency. This is calculated 
using ti MC events, where we evaluate the selection effi- 
ciency separately for twenty bins in lepton 77, to reduce 
the model dependence of our acceptance correction and 
to provide sufficient MC statistics. 

The resolution of the measurement of lepton 77 is ob- 
tained from studies of ti MC events by comparing the 
generated value of 77 with the value measured follow- 
ing event reconstruction. For electrons and muons, we 



use the 77 of tracks measured in the tracking system and 
find this resolution to be the same for both types of lep- 
tons. This resolution is also investigated using cosmic- 
ray muons that appear as dimuon events and is found to 
be w 0.0026, consistent with the MC expectation. For 
« 99.8% of the electrons or muons in ti MC events, the 
sign of lepton 77 is correctly reconstructed. Migration of 
events within the "forward" or "backward" regions does 
not affect the reconstructed angular asymmetry except 
for negligible acceptance corrections. The reconstruc- 
tion effects on the measurement of 77 can therefore be 
neglected for charged leptons. 

The Z+jets background, which is predicted through 
MC simulation contributes to the asymmetry. To 
study the influence of the Z+jets background, we per- 
form measurements of all six asymmetries in a sample 
dominated by Z+jets production in final states with two 
electrons or two muons. Applying the same event selec- 
tions as for the final ti enriched sample, except for the Ifs? 
significance likelihood and IpT requirements, all asymme- 
tries are measured using the same procedure as for the 
measurement of ti asymmetries, but treating Z+jets as 
"signal" and ti as "background" . In this control sample, 
all other background contributions are negligible. The 
data and MC predictions for the 77 distribution of posi- 
tively and negatively charged leptons, for Q ■ 77, and A77, 
are in good agreement, as presented in 21\. 

To verify that the measurement of the ti asymmetries 
is unbiased and correctly estimates the statistical uncer- 
tainty of the result, we perform the measurement using 
ensembles of MC pseudo-experiments. To obtain sam- 
ples with different asymmetries, we mix a ti MC event 
sample weighted to have no asymmetry with different 
fractions of ti MC events with a SM asymmetry. We 
fluctuate the expected number of events in the "forward" 
and "backward" direction for each pseudo-experiment as- 
suming Poisson statistics and apply the same procedure 
as for data to extract the asymmetry. This test shows 
that the measurement is unbiased and that the statisti- 
cal uncertainties are estimated correctly. 

Systematic uncertainties can affect the distributions in 
lepton 77. In particular, the energy scale for jets, jet en- 
ergy resolution, jet reconstruction, the normalization of 
background, the MC-derived acceptance, and the finite 
number of MC events can shift the measured asymme- 
try. The normalization of the background has uncertain- 
ties from diboson and Z+jets cross sections, as well as 
a 6.1% uncertainty on the data sample's integrated lu- 
minosity. The systematic uncertainties on the light and 
heavy-flavor jet energy scales, jet energy resolution, and 
the jet reconstruction can affect the acceptance. We eval- 
uate the size of these uncertainties by applying the varia- 
tion in acceptance corrections and in the differential dis- 
tribution of lepton 77 in deriving the ti asymmetry. 

In addition, we compare the acceptance from single 
leptons obtained from simulated ti events with the accep- 
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TABLE 1: Systematic uncertainties for the six unfolded 
asymmetries defined in Eqs. (|TJ|- (J5j) for the combination of 
all dilepton final states. All values are given in %. 



Source 


A 1 


A l+ 


A 1 " 


A 1 


A u 


A l CP 


Jets 


1.1 


0.8 


1.7 


1.0 


1.5 


1.2 


MC statistics 


0.4 


0.4 


0.4 


0.3 


0.5 


0.3 


Bkg normalization 


0.3 


0.3 


0.6 


0.3 


0.7 


0.3 


Acceptance 


0.7 


0.2 


1.5 


0.7 


2.3 


0.9 


Total 


1.4 


1.1 


2.4 


1.3 


2.9 


1.6 



TABLE 2: Measured asymmetries for leptons, as defined in 
Eqs. H])-©, including statistical and systematic uncertainties 
for the combined dilepton final states using raw and unfolded 
distributions are compared to predictions from mc@nlo in- 
cluding QCD+EW corrections. Our predictions are calcu- 
lated using the NLO QCD+EW distributions in both numer- 
ator and denominator of Eqs. <[TJ) - (J5j . This is different to the 
calculations in Refs. 0, [l7| where the denominator is calcu- 
lated in LO QCD to derive expressions for the asymmetries 
of 0(a a ). All values are given in %. 





Raw 


Unfolded 


Predicted 


A 1 


2.9 ±6.1 ±0.9 


2.5 ±7.1 ±1.4 


4.7 ±0.1 


A l+ 

^FB 


4.5 ± 6.1 ± 1.1 


4.1 ±6.8 ± 1.1 


4.4 ±0.2 


aK 


-1.2 ± 6.1 ± 1.3 


-8.4 ±7.4 ±2.4 


-5.0 ±0.2 


A 1 


3.1 ±4.3 ±0.8 


5.8 ±5.1 ±1.3 


4.7 ±0.1 


A u 


3.3 ± 6.0 ± 1.1 


5.3 ±7.9 ±2.9 


6.2 ±0.2 


A 1 


1.8 ± 4.3 ± 1.0 


-1.8 ±5.1 ±1.6 


-0.3 ±0.1 



tance obtained from Z — > £ + £~ data. We select a data 
sample enriched in Z — > £ + £~ events, where one lepton is 
required to pass tight lepton-selection criteria to function 
as a "tag" and the other "probe" lepton to pass a loose 
lepton selection. The acceptance is evaluated as function 
of rj by applying a tight-lcpton identification requirement 
on the probe. No significant difference is observed be- 
tween the acceptance for positive or negative pseudora- 
pidities, nor between positively and negatively charged 
leptons. A systematic uncertainty on the acceptance is 
defined for each lepton charge by the difference in accep- 
tance between the forward and backward hemisphere of 
the detector. This study is performed separately for elec- 
trons and muons. The systematic uncertainties are added 
in quadrature to yield the total systematic uncertainties 
given in Table [T] 

Using the distributions in Fig. [1] the lepton asym- 
metries of Eqs. (U)-© are measured and corrected for 
acceptance effects ("unfolded"). The measurements of 
the uncorrected (raw) asymmetries as well as the un- 
folded asymmetries are compared to the predictions from 
MC@NLO including QCD±EW corrections @ in Table [2] 
All unfolded asymmetries arc in agreement with the SM 
predictions. 

The asymmetry A FB defined in Eq. ([2]) is also mea- 
sured in ^±jets final states @. The result for A e FB = 
(15.2 ±4.0)% is compared to a predicted value from 



mcOnlo of (2.1 ± 0.1)%. We checked that the predicted 
asymmetry is independent of the final state and the dif- 
ference to our prediction of (4.7 ± 0.1)% is only due to 
the QCD+EW corrections. The dominant systematic un- 
certainty on the prediction and on our measurement in 
dilepton final states is given by jet reconstruction related 
systematics. The total uncertainty of the measurement 
is dominated by the statistical component. Since the 
£±jets and dilepton final states are selected to be statis- 
tically independent, we can improve the uncertainty on 
Ap B by combining both measurements. 

The combination of the two asymmetries Ap B is per- 
formed using the BLUE method 0, [23[. All system- 
atic uncertainties evaluated in both measurements are 
treated as fully correlated. The combination yields a 
leptonic FB asymmetry of A FB = (11.8 ± 3.2)%, where 
the £±jets channel contributes 63.9% and the dilepton 
channel 36.1% of the information. This represents an 
improvement of about 20% relative to the uncertainty 
in the £+jets channel alone. The consistency between 
the two individual measurements is 68%. Comparing the 
combined result to the predicted leptonic FB asymmetry 
from MC@NLO plus higher order QCD+EW corrections, 
A FB (predicted) = (4.7 ± 0.1)%, we observe a disagree- 
ment at the level of 2.2 SD. 

To further investigate this deviation of the asymmetry 
from the SM prediction, we analyze the longitudinal po- 
larization of the top quark. While in the SM top quarks 
are expected to be produced unpolarized in tt events, 
there are many beyond the SM models that would en- 
hance the tt FB asymmetry [l| and therefore the lep- 
tonic asymmetries defined in Eqs. H])-©, and would 
also lead to a non-vanishing longitudinal polarization of 
the top quark. Examples are models with new parity- 
violating interactions. In the absence of effects from ac- 
ceptance, the distribution of cos8~ and cos# + should be 
isotropic [|| for unpolarized top quarks, where 9 + (9~) 
is the angle between the direction of the £ + (£~ ) in the t 
(t) rest frame and the t (t) direction in the tt rest frame. 
A longitudinal polarization of the top quark would cause 
asymmetric cos^* distributions. 

Assuming CP invariance, i.e. that the distributions 
of cos 9 + and cos 9~ are equal, we measure the distribu- 
tion cos 9, defined by the sum of the cos ± distributions. 
The calculation of the angles 9 requires a transforma- 
tion of the momenta of the charged leptons into the t 
and t quark rest frames. Every event must therefore be 
fully reconstructed. This is performed using the neutrino 
weighting method, devised originally to measure the top 
quark mass in the dilepton channel 12 41 and recently ap- 
plied to measure tt spin correlations [201 ] - 

In Fig. [2j the cos 9 distribution is shown separately for 
the dilepton and ^+jets final states. The distribution for 
tt events produced via a leptophobic topcolor Z' boson, 
with the same parity- violating couplings to quarks as the 
SM Z boson and a width T = 0.012M Z ml [51 is also 
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FIG. 2: The distribution of cos 6 is shown for the combination 
of the dilepton channels (a) and the ^+jets channels (b). The 
data are compared to the SM predictions. The vertical error- 
bars on the data points indicate the statistical uncertainty of 
the data. The distribution of ti pairs produced via a hypo- 
thetical Z' boson is also shown; the uncertainty due to the 
limited size of the MC sample is shown by the shaded band. 
The same Z' model as in [29, ^al is used. 



shown to illustrate the effect of producing t and t quarks 
with longitudinal polarization. The agreement between 
the data and the SM prediction in both distributions is 
good, yielding a Kolmogorov-Smirnov test probability of 
14% in the dilepton channel and 58% in the ^+jets chan- 
nel. There is no significant hint of new sources of parity 
violation leading to a longitudinal polarization in ti pro- 
duction. 

In conclusion, we measured angular asymmetries in ti 
production based on rj distributions of charged leptons. 
We find the leptonic FB asymmetry Ap B and the lepton 
asymmetry A u in agreement with the SM prediction in 
the dilepton final state. Combining our measurement of 
Ap B with the measurement performed using leptons in 
£+jets final states yields Ap B = (11.8 ± 3.2)%, which is 
2.2 SD above the higher order QCD+EW prediction of 
Af B (predicted) = (4.7 ± 0.1)%. The top-quark polar- 
ization in the dilepton and £+jets final states show good 
agreement between the data and the SM prediction. 
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and negatively (b) charged leptons, the distribution of q ■ r\ (c) and the distribution of A?7 = rj e + — r\ t - (d) are shown. 
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